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Abstract 

The  discontinuous  Galerkin  (DG)  method  continues  to  maintain  heightened  levels  of  interest 
within  the  simulation  community  because  of  the  discretization  flexibility  it  provides. 
Although  one  of  the  fundamental  properties  of  the  DG  methodology  and  arguably  its  most 
powerful  property  is  the  ability  to  combine  high-order  discretizations  on  an  inter-element 
level  while  allowing  discontinuities  between  elements,  this  flexibility  generates  a  plethora  of 
difficulties.  First,  by  evolving  more  degrees  of  freedom,  a  more  restrictive  cfl  condition  is 
needed,  which  is  computationally  expensive.  Secondly,  the  lack  of  continuity  across  element 
interfaces  hampers  visualization  efforts  as  many  visualization  packages  assume  higher  levels 
of  continuity.  However,  DG  also  has  the  property  of  superconvergence.  That  is,  it  achieves 
essentially  twice  the  usual  convergence  rate  at  specific  points  within  the  mesh,  or  in  a  special 
norm,  allowing  for  faster  convergence  to  a  reasonable  solution. 

The  goal  of  this  research  is  to  exploit  the  inherent  property  of  superconvergence  in  order  to 
(1)  create  a  better  pairing  between  DG  solutions  and  the  time-stepping  method  used,  and  (2) 
create  filters  that  allow  highly  accurate  visualization  of  DG  data. 

Portions  of  this  work  was  performed  in  active  collaboration  with  Prof.  Robert  M.  Kirby  at  the 
University  of  Utah,  who  is  sponsored  by  the  Air  Force  Office  of  Scientific  Research,  Air 
Force  Material  Command,  USAF,  under  grant  number  FA9550-12-1-0428. 

Status/Progress 

AFOSR  funding  to  support  this  research  was  obtained  in  February  2013.  This  funding  was 
used  to  support  Mr.  Xiaozhou  Li  (TU  Delft,  February  2013  -  July  2015),  Mr.  Daniel  Frean 
(UEA,  July  2013  -  July  2016)  and  Ms.  Julia  Docampo  (October  2013  -  July  2016).  Mr.  Li 
successfully  defended  his  thesis  in  July  2015.  Mr.  Frean  and  Ms.  Docampo  are  expected  to 
submit  their  thesis  in  October  of  2016.  Additionally,  this  funding  was  used  to  support 
collaborative  visits  related  to  this  grant.  Ms.  Docampo  visited  the  University  of  Utah  to  work 
with  Prof.  Robert  M.  Kirby  in  July  of  2014  and  February-March  of  2016.  Mr.  Frean  visited 
the  University  of  Massachusetts  -  Dartmouth  to  work  with  Prof.  Sigal  Gottlieb  in  February- 
March  of 2016. 

This  funding  has  allowed  the  PI  to  make  several  contributions  over  the  lifetime  of  this  grant: 

•  In  the  first  year  of  this  grant,  the  main  accomplishments  were:  1)  demonstrating 
numerically  that  there  is  indeed  a  “sweet  spot”  for  scaling  of  the  Smoothness- 
Increasing  Accuracy  Conserving  (SIAC)  filter  when  applied  to  nonunifonn  meshes; 

2)  Theoretically  and  numerical  demonstration  of  the  2k+l  order  accuracy  of  the  SIAC 
filter  when  applied  to  meshes  with  smoothly-varying  mesh  sizes;  3)  Numerically 
demonstrating  the  usefulness  of  the  boundary  SIAC  filter  for  nonlinear  equations  [4- 
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7,9,15-17,20];  and  4)  Linking  the  known  results  for  the  eigenstructure  analysis  of  the 
discontinuous  Galerkin  (DG)  method  for  the  upwind  flux. 

•  In  the  second  year  of  the  grant,  the  main  accomplishments  were:  1)  Establishing  a 
more  theoretical  and  numerical  understanding  of  a  computationally  efficient  scaling 
for  the  SIAC  filter  for  nonuniform  meshes  [7];  2)  Establishing  connections  to  signal 
processing  and  generalizing  the  SIAC  kernel  through  quasi-interpolation  [19];  3) 
Introducing  a  Hexagon  SIAC  filter  to  reduce  the  kernel  footprint;  4)  Establishing  the 
theoretical  viability  of  the  SIAC  filter  for  nonlinear  scalar  hyperbolic  conservation 
laws,  which  includes  examining  the  role  of  divided  difference  estimates  [18,25-26];  5) 
Establishing  superconvergence  of  Semi-Lagrangian  DG  (SLDG)  methods  [12];  as 
well  as  6)  Establishing  theoretical  and  numerically  that  the  upwind  biased  (non¬ 
monotonic)  flux  is  more  beneficial  for  even  polynomial  degree  approximations  than 
the  upwind  flux.  This  was  used  to  computationally  demonstrate  the  role  of  the  flux  in 
minimizing  errors  in  long-time  simulations  [2,23]. 

•  In  the  last  year  of  the  grant,  there  have  been  two  major  areas  of  research:  1)  Fully 
discrete  estimates  for  the  dispersion  and  dissipation  errors  were  obtained  for  various 
time-stepping  strategies  [3,24];  and  2)  Establishing  the  optimal  kernel  footprint  for 
the  SIAC  filter  for  computational  efficiency  [1,13,22].  The  key  findings  for  each  of 
these  areas  are  summarized  in  the  next  section. 

Key  Findings  in  final  year  of  grant 

•  Fully  discrete  estimates  for  the  dispersion  and  dissipation  errors.  For  this  study  the 
following  time-stepping  schemes  were  paired  with  the  DG  scheme:  the  Runge-Kutta 
scheme  of  2nd,  3ld  and  4th  orders  (RK2,  RK3,  RK4),  Lax-Wendroff  (LW),  New  Lax- 
Wendroff  (nLW)  and  Two  Derivative  Runge-Kutta  (TDRK)  schemes.  For  LW  and 
TDRK,  the  dissipation  and  dispersion  errors  are  only  k+1  order  accurate.  For  the  RK 
schemes  and  nLW  scheme,  the  dispersion  error  is  superconvergent  of  2k+l.  The 
nLW  scheme  imposes  an  auxiliary  formulation  for  the  derivatives  in  order  to  recover 
the  2k+l  order  accurate  dispersion  error.  This  technique  was  introduced  for  the 
TDRK  (labelled  nTDRK)  and  the  superconvergent  dispersion  error  was  obtained. 

The  equivalence  in  long  time  error  behaviour  for  the  RK  schemes  as  well  as  nLW  and 
nTDRK  was  established  (Figure  1)  [?]. 


Figure  1:  Time  history  plots  of  the  DG  spatial  discretization  paired  with  the  Runge-Kutta  (RK),  Lax- 
Wendroff  (LW)  and  Two  Derivative  Runge-Kutta  (TDRK)  schemes  for  k=l  (Left  two  plots)  and  k=2 
(Right  two  plots).  The  solid  blue  and  black  lines  represent  the  original  LW  and  TDRK  schemes.  The 
dashed  red  line  represents  and  upwind-biased  flux  within  the  DG  scheme.  For  k=l,  the  upwind  DG  with 
RK  is  equivalent  to  nLW  and  nTDRK.  For  k=2 ,  the  upwind-biased  flux  paired  with  the  Runge-Kutta 
scheme  is  optimal. 

•  Establishing  the  optimal  kernel  footprint  for  the  SIAC  filter  for  computational 

efficiency.  This  endeavour  consisted  of  two  key  components:  1)  the  effect  of  rotating 
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the  kernel;  and  2)  numerical  and  theoretical  demonstration  of  the  ability  to  use  a  one¬ 
dimensional  SIAC  kernel  for  multi-dimensional  data.  In  computational  experiments, 
it  was  noted  that  given  the  correct  kernel  scaling  for  the  rotated  kernel, 
superconvergence  in  the  errors  as  well  as  reduction  in  the  errors  occurred.  However, 
the  computational  efficiency  of  rotating  the  kernel  hampered  simulations  because  of 
the  necessity  to  determine  the  intersections  between  the  DG  mesh  and  kernel  mesh. 
When  the  theory  was  investigated,  it  was  detennined  that  a  tensor  product  filter  was 
not  necessary  for  multi-dimensional  data.  The  error  between  the  exact  solution  and 
SIAC  filtered  solution  remains 

||  u-Uh  ||  <  =  Ch2p+1. 

This  presents  a  significant  reduction  in  the  kernel  footprint  as  well  as  an  increase  in 
computational  efficiency  (Figure  2). 


z 

E 

; ; 

Figure  2:  Illustrations  of  the  DG  mesh/Kernel  mesh  intersection  regions.  Each  region  requires 
implementation  of  a  multi-dimensional  quadrature  rule.  Left:  Cartesian-coordinate  aligned  filter.  Center: 
Rotated  filter.  Right:  One-dimensional  filter. 

The  one-dimensional  kernel  for  two  rotations,  p/4  and  3p/4,  was  then  tested  on  the 
DG  approximation  to  a  two-dimensional  advection  equation.  This  was  compared 
with  the  traditional  cartesian-aligned  tensor  product  filter.  Results  demonstrated  error 
reduction  and  superconvergence  (Table  1).  That  is,  the  SIAC  Line  filter  is  able  to 
improve  the  convergence  rate  of  the  DG  approximation  from  k+1  to  2k+l  while 
reducing  the  error.  The  contour  plots  for  the  errors  are  given  in  Figure  3. 
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L2-e  rror  Order 

L2-e  rror 
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20 
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- 
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- 
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- 

1.6e-03 

40 

2.4e-03 

2.02 
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2.98 

2.6e-04 

3.33 

2.0e-04 

80 

5.9e-04 

2.01 

2.4e-05 

2.99 

2.8e-05 

3.21 

P2 

20 

2.4e-04 

- 

1.5e-06 

- 

1.4e-04 

- 

6.1e-06 

40 

2.9e-05 

3.01 

4.7e-08 

4.99 

2.3e-06 

5.91 

1.2e-07 

80 

3.6e-06 

3.01 

1.5e-09 

5.00 

3.7e-08 

5.95 

- 

p3 

20 

4.5e-06 

- 

7.7e-10 

- 

1.6e-05 

- 

1.4e-07 

40 

2.8e-07 

4.01 

6.9e-12 

6.79 

6.9e-08 

7.87 

5.6e-10 

Table  1:  Results  of  the  one-dimensional  SIAC  filter  applied  to  the  DG  approximation  to  a  2D  linear 
advection  equation  with  sine  initial  condition.  The  left  two  columns  represent  the  errors  and  order  in 
the  DG  solution,  the  right  column  the  error  of  the  traditional  tensor  product  SIAC  filter  and  the  middle 
four  columns  the  errors  and  orders  of  the  one-dimensional  SIAC  filter  for  two  different  rotations.  The 
one-dimensional  SIAC  filter  reduces  the  error  and  increases  the  order  of  accuracy  from  k+1  to  2k+l. 
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Figure  3:  Contour  plots  of  the  error  for  the  DG  approximation  (left),  tensor  product  SIAC  filter 
(center  left),  SIAC  line  filter  with  a  rotation  of  p/4  (center  right),  and  SIAC  line  filter  with  a  rotation  of 
3p/4  (right)  for  the  2D  advection  equation  with  sine  initial  condition. 


Most  Significant  Accomplishments 

Over  the  period  of  the  grant,  the  most  significant  accomplishments  are: 

(1)  Establishment  of  the  relationship  of  the  flux  in  the  discontinuous  Galerkin 
approximation  to  superconvergent  dissipation  and  dispersion  errors.  This  ensures 
that  waves  are  correctly  propagated  forward  in  time,  and  thus  ensuring  a  reduced 
long-time  error  [23-24]. 

(2)  Establishing  the  SIAC  line  filter  both  computationally  and  theoretically.  Using  a 
one-dimensional  filter  for  multi-dimensional  data  has  a  reduced  computational 
cost  that  ensures  error  reduction  as  well  as  an  improvement  in  the  convergence 
rate.  Additionally  it  ensures  smoothness  of  the  approximation  in  all  directions 
[22]. 

(3)  Establishing  connections  to  signal  processing  by  investigating  the  relation 
between  SIAC  filtering  and  quasi-interpolation  [19]. 

(4)  Establishment  of  a  computationally  efficient  boundary  SIAC  filter.  This  allows 
for  an  increase  of  the  smoothness  of  the  DG  approximation  and  generally  a 
reduction  in  the  error  [15-17,20]. 
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